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ABSTRACT 

We report new imaging polarimetry observations of the Galactic compact HII region K3-50 
using CanariCam at the Gran Telescopio Canarias. We use a standard polarimetric analysis 
technique, first outlined by Aitken, to decompose the observed polarisation images centred 
at 8.7, 10.3, and 12.5 41 m into the emissive and absorptive components from silicate grains 
that are aligned with the local magnetic field. These components reveal the spatially-resolved 
magnetic field structures across the mid-infrared emission area of K3-50. We examine these 
structures and show that they are consistent with previously observed features and physical 
models of K3-50, such as the molecular torus and the ionised outflow. We propose a 3D ge¬ 
ometry for all the structures seen at different wavelengths. We also compute relevant physical 
quantities in order to estimate the associated magnetic field strengths that would be implied 
under various physical assumptions. We compare these results with MHD simulations of pro¬ 
tostar formation that predict the magnetic field strength and configuration. We find that the 
magnetic field may be dynamically important in the innermost 0.2 pc of the molecular torus, 
but that the torus is more likely to be rotationally-supported against gravity outside this ra¬ 
dius. Similarly, magnetic fields are unlikely to dominate the global physics of the ionised 
outflow, but they may be important in helping confine the flow near the cavity wall in some 
locations. Ours is the first application of the Aitken technique to spatially-resolved magnetic 
field structures in multiple layers along the line of sight, effectively a method of “polarisation 
tomography.” 

Key words: stars: circumstellar matter — techniques: polarimetric — ISM: magnetic fields 
— ISM: individual (K3-50) — stars: formation. 


1 INTRODUCTION 

Magnetic fields may play a number of important roles in various as- 
trophysical settings, such as in AGN (synchrotron emission, plasma 
confinement), galactic spiral arms (large-scale instabilities that col¬ 
lect material into molecular clouds), or the formation of stars or 
planets (pressure support against gravity, angular momentum loss). 
However, whether magnetic fields have a dominant role in these 
settings, especially during star formation, is still unclear, and so the 
study of magnetic fields in such objects is a topic of active research. 
Observationally, magnetic fields are difficult to measure, 
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and no single technique can simultaneously give the full three- 
dimensional field strength and orientation. However, significant 
partial information can be obtained with certain techniques. Polari¬ 
metric photometry gives the magnetic field orientation projected 
onto the plane of the sky, but without direct information on the 
field strength or the line-of-sight contribution (e.g., |Aitken |1996| 
|Smith et al.|2000{ |Aitken et al.|2 004). In contrast, Zeeman spec¬ 
tral line measurements give the line-of-sight field strength (a one¬ 
dimensional projection), but are much less sensitive to the orthog¬ 
onal components’ strength or orientation (see review by |Crutcher| 

[20l2| . 

Both techniques are technically challenging, since in both 
cases magnetic field measurements depend on obtaining accurate 
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values for the polarised contributions to emission or absorption 
(e.g., the Stokes parameters Q, U, and/or V ), which are usually 
quite small compared to the total intensity I (although scattering in 
the near-infrared can give polarisations up to 100%). In the case of 
polarimetry, there is the additional complication of identifying the 
source of the polarisation, which (depending on wavelength and 
particle size) can be due to Rayleigh, Mie, or electron scattering, or 
dichroic absorption and/or emission from dust grains. 

Thus, experiments capable of yielding high quality, spatially 
resolved images of either magnetic component (i.e., plane of sky 
or line of sight) are highly prized, not only because they are still 
rare, but also since some measure of the three-dimensional mag¬ 
netic field structure can be recovered from such data by context. 

As part of a multi-source campaign to better understand the 
role of magnetic fields in astrophysical disks in all settings, we 
report here one of the first results from the newly-commissioned 
mid-infrared (MIR) CanariCam instrument, a multi-mode user fa¬ 
cility for the Gran Telescopio Canarias (GTC). We used CanariCam 
in its imaging polarimetry mode with a series of MIR filters to ob¬ 
serve the polarisation of the 10 /rm silicate feature towards the mas¬ 
sive star-forming region K3-50, which |Wright| ( j2007) selected as 
an ideal candidate for high-resolution MIR polarimetry. Our intent 
was to use these observations to construct, from the dichroic prop¬ 
erties of silicate grains in this source, one of the first examples of 
a spatially-resolved map of the magnetic field structure, from both 
emitting and absorbing layers along the line of sight and across 
the field of view. We wanted to use this “polarisation tomography” 
to elucidate the nature and configuration of this source’s magnetic 
fields, and determine whether any detected fields in the disk, out¬ 
flow, or embedded objects were playing a dominant role (or not) in 
the evolution of this region. 

K3-50 is part of the radio source W58, located in the outer 
Galaxy at (l,b) = (70°3,+l °6). The distance to K3-50 has been taken 
as 8.7 kpc in most recent work, as originally determined by Harris 
< 11975) . However, |Palmer et al.| ( |2003l > pointed out that this should 
be rescaled to 7.4kpc with a value of Ro = 8.5 kpc for the Sun’s 
Galactocentric distance, and proportionate rescalings would apply 
for other values (e.g., d = 7.3 kpc with a more modem Ro = 8.4 kpc; 
|Reid et al.|2009). While acknowledging this, we nonetheless adopt 
d = 8.7 kpc to be consistent with Ho ward et al.| <fl997 1 in particular, 
so that we can compare our results with their model for the molec¬ 
ular toms ( §4.2. 1) . A 16% smaller distance will obviously reduce 
all masses and luminosities discussed herein by 29%. 

K3-50 is an archetypical luminous (~2x 10 b L©) compact HII 
region with ongoing massive star formation, and has been observed 
by many workers at multiple wavelengths. For the purposes of this 
study, relevant prior results include the 2 cm continuum and recom¬ 
bination line interferometric study of the larger area around K3- 
50 by |De Pree et al.[ < jl994) , showing the existence of a narrow, 
ionised, bipolar outflow emanating from the central ultracompact 
HII region (UCHII) K3-50A; the 3 mm molecular interferometry 
maps of |Howard et al ,| ( j 1 997) , showing an inclined toroid or disk 
of size ~0.4pc with a central cavity coincident with the K3-50A 
UCHII peak; and the near-IR speckle interferometry of jHofmann| 
[et al.||2004) , showing a number of bright knots inside a biconical 
nebulosity at the centre of the molecular toroid’s cavity. In addition, 
|Okamoto et al.| ( |2003[ > obtained MIR imaging and spectroscopy of 
the HII region’s ionising sources in the central cluster. 

Because of these geometrical features and its brightness at 
MIR wavelengths, K3-50 makes an ideal object for polarisation 
mapping of its magnetic field configuration in the 10 gm silicate 
feature. Besides the new data and analysis presented here, the only 


other published IR polarimetry of K3-50 IRS were presented by 
|Dyck & Capps|j 1 978) and|Smith et al7|(|2000| >, with a non-detection 
at 100/rm by |Novak et al.|(l989) . |Smithet al.| ( [2000] l also performed 
a decomposition of the absorptive and emissive components at their 
5” resolution, and an analysis of the orientation of these compo¬ 
nents in terms of disk/toroid and outflow axes was published by 
|Wright[ < |2007) . 

We describe the camera and observational details in the next 
section, and present the results of the data reduction and analysis in 
§3. We discuss the implications of these results in §4, and end with 
a brief conclusion. 

2 OBSERVATIONS 

2.1 CanariCam 

CanariCam is a multi-mode instrument for the GTC, offering imag¬ 
ing and spectroscopy in the mid-infrared N (7.5-13.5 (im) and Q 
(16-26 /im) atmospheric bands, and polarimetry in the N band. De¬ 
scriptions of the overall instrument design and of the polarimetric 
mode have been given by Telesco et al. ( 2003) and Packha m et al.| 
( 2005) ; see also http://www.gtc.iac.es/instruments/canaricam/cana- 
ricam.php. CanariCam is built around a 320x240 Raytheon array, 
forming images at a pixel scale of 0'/08 and giving a native field of 
view of 25”6x 19” 2. CanariCam is equipped with a range of broad- 
and narrow-band filters with central wavelengths between 7.8 and 
24.5 //m. The theoretical angular resolution of the camera at 10 /rm 
is 0725 (1.22A /D), but the seeing at the GTC degrades this some¬ 
what, depending on observing conditions (see §2.2) . 

During integration in the polarimetry mode, the light beam 
passes though a CdSe half-wave-plate (HWP) modulator, which is 
rotated sequentially to four orientations (0°, 22?5, 45°, and 67?5), 
to allow for full polarisation calibration during data reduction. 
Then, a CdSe Wollaston prism is used to separate the main opti¬ 
cal beam into two orthogonally polarised ones (i.e., the o ray and 
e ray). At the stop of each HWP orientation, a frame containing 
both o ray and e ray images is recorded simultaneously. Depend¬ 
ing on the integration time, several tens to hundreds of frames are 
obtained, and saved in a multi-dimensional FITS file at the end of 
each integration. 

This dual-beam configuration minimises the effect of seeing 
and “common-mode noise,” such as that from atmospheric vari¬ 
ations, and is the first ever dual-beam MIR polarimeter, offering 
significant gains in sensitivity over previous instruments (over and 
above being on a 10 m telescope). To prevent overlap between the 
o ray and e ray images of an extended object, a focal plane mask is 
used, and the effective field of view reduces to three narrow rectan¬ 
gular apertures, each of which is approximately 25 , /6x3 , /0. 

2.2 K3-50 Data Collection, Calibration, and Reduction 

We observed K3-50 on UT 2012 July 6 and 8 as part of the Canari¬ 
Cam Science Team guaranteed time for the GTC. We positioned 
the CanariCam field of view over the brightest part of the near- 
IR speckle image of |Hofmann et al.| < |2004) , oriented with the de¬ 
tector’s long axis exactly N-S. The three filters Si2, Si4, and Si6 
were used, centred at 8.7, 10.3, and 12.5 /rm with bandwidths 1.1, 
0.9, and 0.7 /j,m, respectively. This combination of filters allows us 
to separate the dichroic absorption and emission components from 
the total MIR polarisation using the methodology of |Aitken et al.[ 
(20041 (see §3.1|below). The observations of K3-50 were interlaced 
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Figure 1. Colour composite image of the 8.7 (im (Si2, blue), 10.3 //m (Si4, 
green), and 12.5 (im (Si6, red) I images. The brightness scales in each I 
image have been adjusted to maximise the colour contrast in this image, 
which then represent the relative variation in flux ratios between the 3 fil¬ 
ters. Thus, the eastern peak has relatively higher Si6/Si4 and Si6/Si2 ratios 
than other parts of the image, and so appears redder; the western peak is 
relatively brighter in Si2 & Si4 compared to the rest of the image, and so 
appears as a cyan colour; the broader area between and around these peaks 
has a greener hue, indicating relatively high Si4/Si2 and Si4/Si6 ratios in 
this area; and the far western edge of the western peak is almost pure blue, 
indicating relatively high Si2/Si4 and Si2/Si6 ratios compared to other ar¬ 
eas. These colours/ratios persist to fainter flux levels as well: red ratios to 
the east, cyan to the north, blue to the west, and green to the south includ¬ 
ing the “hook” feature. Contours in percent of the peak of the Si4 image are 
shown at 0.8,1.5, 2.6, 3.8, 6, 12 (pale yellow), and 25 to 95 by 10 (magenta). 


with 100 s on-source integrations on a nearby PSF calibrator HD 
168723 and a polarisation/PA calibrator AFGL 2591. The former 
is a Cohen standard star (| Cohen et al. | |1999 ), while the latter is 
chosen from the atlas of |Smith et al.[j2000fr . The effective reso¬ 
lution measured on the PSF star was found to be ~0”3, indicat¬ 


ing a moderate influence from seeing, but the image quality is still 
diffraction-dominated. 

For K3-50, the total on-source integration time was 864, 437, 
and 728 seconds (resp.) for each filter. The resulting rms noise in 
the (unfiltered, unsmoothed, unbinned, raw) total intensity I im¬ 
ages is 13, 21, and 21 mjy arcsec -2 . We used the standard chop- 
nod technique to correct for the thermal background from the sky 
and the telescope. The chop throw was 20” along the E-W direc¬ 
tion and the chop frequency was 2 Hz. No flat-field correction is 
needed for instruments like CanariCam, since the pixel response is 
very uniform across the entire array. 

We present a 3-colour composite of the raw I images in Fig¬ 
ure [I] where one can see by the rendered colours that different 
structures in K3-50 have different flux ratios among the 3 filters. In 
this rendering, redder colours probably indicate cooler dust and/or 
higher extinction, while bluer colours may indicate warmer and/or 
less obscured features. 

To recover the Stokes parameters from the raw data, an IDL- 
based custom software package iDealCam jLi et al.||20 13} was 
used. The procedures are briefly summarised here. After first-order 
background correction by chopping, a pair of sub-images with 
identical dimensions (6”6x2"8) was extracted from each frame, 
corresponding to the o ray image and the e ray image of the source. 
The sum of the two sub-images gives a measure of Stokes / (total 
intensity), and the difference between them is a measure of Stokes 
Q or U, depending on the HWP angle. All frames in the raw FITS 
file are processed in the same way, and the results are averaged 
at the end to generate the final images of /, Q , and U. The stan¬ 
dard deviation between frames is also computed to estimate the 
uncertainty of each Stokes parameter. For the purpose of data visu¬ 
alisation, we also calculated the (linear) polarisation percentage P 
(defined as \JQd + U 2 / 1) and the angle of polarisation 9 (defined 
as 0.5 tan' 1 [U/Q]) for each pixel (or binned pixel to enhance the 
signal-to-noise ratio) and at each wavelength. Finally, the PA cali¬ 
brator was used to translate 9 to the reference system with respect 
to the true north on the sky. Typical uncertainties in these quantities 
for a 2x2-binned pixel range from ~0.1 to 0.5% in P, and ~5-15° 
in 9 , for the brightest and outlying areas of emission, respectively. 
These products are computed in iDealCam by formal propagation 
of measurement uncertainties in /, Q, and U, and are shown in 
combined form in Figure [2] 

Since CanariCam is installed at the Nasmyth focus of the tele¬ 
scope, the tertiary flat mirror introduces significant (~0.6%) in¬ 
strumental polarisation (P-inst). This instrumental effect has been 
characterised during the commissioning of CanariCam, and its be¬ 
haviour is well known, so P-inst was corrected during data reduc¬ 
tion, including allowing for the different parallactic angle between 
the source and PA calibrator. A detailed description of the data col¬ 
lection and data reduction procedures, including P-inst correction, 
will be found in Li et al. 2015 (in prep.). 


3 DATA ANALYSIS 

3.1 Polarisation Components — The Aitken Method 

Aitken and co-workers have shown that the 10 fim spectral sig¬ 
nature of silicate particles provides a powerful diagnostic tool for 
probing the complex polarimetric properties of astrophysical envi¬ 
ronments where silicates are a significant component ( jAitken|1996| 
|Smith et al.|2000||Aitken et al.|2004| >. This is evident in Figure |2J 
where one can see a much stronger variation from panel to panel 
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Figure 2. CanariCam imaging products for K3-50 at (left to right) 8.7, 10.3. and 12.5 fi m, derived from the corresponding I, Q, U images using iDealCam. 
Note that these panels are truncated in their long dimension and a S/N > 3 clip applied to show only the high significance data. In each panel, the greyscale 
image is fo r the total intensity I, in units indicated by the corresponding greyscale wedge at the top; these are consistent with images presented by |Okamoto| 
[et al. (2003} and will not be further discussed here. The blue contours are for the polarised intensity (Si2, first at 150mJy arcsec -2 , rising in steps of 
75 mJy arcsec -2 ; Si4, 300(150)mJy arcsec -2 ; Si6, 300(120) mJy arcsec -2 ); and the red vectors show the polarisation orientation and percentage as indicated 
by the red scale bar at the bottom. The FWHM of the camera+telescope PSF (including seeing) is shown for reference as a hatched circle in each panel; these 
have sizes 0727, 0731, and 0735, resp. These images have been formed with a 2x2 binning from the original data, so the pixel scale here is 0716. 


in the polarisation maps (blue contours and red vectors) than in the 
greyscale I images. In essence, Aitken et al. showed that, for small 
polarisations, the observed polarisation as a function of wavelength 
can be approximated by a linear combination of the polarisation 
due to dichroic absorption and that due to dichroic emission. They 
then assumed that the spectral profile / of the grain polarisation 
obeys / em ( A) ~ / a bs(A)/r(A), and that the measured absorption 
profile t( A) for the BN object in Orion is typical of silicate grains 
in all sources. With these assumptions, the absorptive and emissive 
components reach their maximum values at 10.3 /im and 11.5 /im, 
respectively. 

This method does not assess the contribution of any other po¬ 
larisation processes, like MIR dust or electron scattering, to the 
data. Instead, the working assumption that such contributions are 
small compared to these dichroic effects, is verified (or not) by the 
goodness of fit of the decomposition to the polarisation data. 


With this approach, all that is needed to solve for the separate 
absorption and emission contributions is a full set of polarisation 
measurements for at least two wavelengths across the silicate ab¬ 
sorption feature, as explained by |Aitken et al.||2Q04[ l. The fact that 
the solutions for P a b s and P em often successfully reproduce the ob¬ 
served polarisations in a variety of sources (e.g. |Smith et al.|2000[ > 
underscores the reasonableness of the assumptions built into this 
technique (including the use of the BN profile), even when multi¬ 
band polarimetry (or spectropolarimetry) is available (i.e., when the 
solutions are over-determined). 


Using the Aitken method in iDealCam, with a functional form 
for t(A) given, e.g., by Gillett et al. (1975) for pure silicate grains, 
we derived from the data shown in Figure[2]the corresponding solu¬ 
tions for the absorptive P a b s and emissive P em components of the 
polarisation, after first convolving the shorter-wavelength images 
to the 12.5/rm resolution. These solutions are shown in Figure[3]in 
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Figure 3. Aitken method decomposition of polarisation components for K3-50, now computed for each 3x3 binned pixel to improve the S/N in this polari- 
metric decomposition. The absorptive component is shown in the left panel, the emissive on the right. In each panel, the red vectors show the polarisation 
orientation and percentage for the respective component; the green vectors within the red vectors indicate the relative uncertainty in the polarisation per¬ 
centages, but are capped to be no longer than the red vectors. The two black-and-white dotted circles in each panel show the locations of the two sample 
decompositions in Fig.[4] The absorptive polarisation reaches a maximum of 10% at a positional offset (Aa,AS) near (-0 / /l,-0 / /2), while the emissive polar¬ 
isation also peaks at 10% near (-0 / /2,-l / /8), counting only higher S/N points within the lowest blue contour. The greyscale image, blue contours, and other 
details are the same in both panels, and the same as in the 10.3 /im panel in Fig.[2] At a distance of 8.7 kpc, the physical scale is = 8700 AU = 0.042pc. 
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Figure 4. Sample spectral decompositions on polarimetric measurements at three wavelengths, using the Aitken method. The two spots chosen for this 
illustration are at the positions indicated in Fig.[3] one near the central region of K3-50, and the other on the outflow cavity wall. Similar to the presentation in 
| Aitken et al. 1(2004) , decomposition results are given in spectra of total polarization (p ), position angle ( 0 ), Stokes Q/I and U/I. The solid line shows the best 
fit from a combination of absorptive and emissive components, which are plotted as dotted and dashed lines, respectively. 
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similar format to Figure [2] with graphical examples of the decom¬ 
position shown in FigureBlfor the two circled positions indicated in 
Figure [3] Note that these solutions show large, coherent structures 
across the field of view, over a scale of many independent resolu¬ 
tion elements. Such features should be indicative of the inherent 
properties of the source, given the comparatively low measurement 
uncertainties. 

Similarly to |Smith et al.H2000) and |Aitken et al ,| ( f2004[ i, we 
also experimented with fitting single (emissive or absorptive) com¬ 
ponents to the CanariCam data (as opposed to fitting both com¬ 
ponents, as described above). Using the \ 2 value of the fit as a 
general criterion complemented by visual examination, we found 
that single-component fits are favoured for only a small portion of 
pixels (<10%). For such pixels, replacing the results from two- 
component fitting with single-component results causes negligible 
effects on the spatial structures observed, in either the emissive or 
absorptive polarisation maps. Therefore, for clarity and simplicity, 
only the two-component fitting results are presented here. 


3.2 Magnetic Field Orientation and Polarisation 
Tomography 

The Aitken decomposition is significant because the absorption and 
emission of the polarised radiation is tied to the physical orien¬ 
tation of silicate grains aligned in an embedding magnetic field. 
For the absorption, the polarisation direction should be the same as 
the magnetic field orientation (projected onto the plane of the sky). 
Since the absorption is greater in the direction parallel to the grains’ 
long axis, which is expected to have a net orientation perpendicu¬ 
lar to the magnetic field, the absorptive polarisation .E-vector will 
be parallel to the field direction. In contrast, the polarised emission 
that we see is preferentially oriented perpendicular to the magnetic 
field ( |Lazarian|2007| [Lazarian & Floang 2014). Therefore, to ex¬ 
amine the intrinsic magnetic field directions in the plane of the sky 
in the absorbing and emitting layers, we can overlay the two pan¬ 
els of Figure [5] but with the emissive vectors rotated by 90°. This 
is shown in Figure [5] where a very interesting trend emerges on 
the relative orientation of the two magnetic field components. To 
quantify and further explore this trend, we plot the position angle 
difference between the emissive and absorptive magnetic field lay¬ 
ers in Figure[6] 

From Figure[6]we see that there are many pixels (filled circles) 
where the magnetic field orientation from both components is well- 
aligned (A 6b < 20°). If there are widespread areas where both 
the emission and absorption from silicate grains seem to trace the 
same magnetic field, it is possible that this field is uniform across 
absorption and emission regions separated along the line of sight, 
or that these polarisation signatures are collocated or commingled, 
perhaps in a clumpy medium threaded by a common magnetic field. 

In contrast, for areas where the polarisation indicates differ¬ 
ently aligned magnetic fields (A 9b > 20°, open circles), it is rea¬ 
sonable to suppose that the two magnetic field components are lo¬ 
cated in two regions separated along the line of sight. In this case, 
we further suppose that the absorption-related magnetic field lies 
in front of the emission-related one; it is then also likely that the 
emission-related field, being associated with warm dust, lies closer 
to the main heating source for the region. 

The distribution of points in Figure [6] show some clear order¬ 
ing, and suggests that the pixels with similar magnetic field char¬ 
acteristics may also be spatially correlated with each other. To ex¬ 
amine this, we construct three subsets of the magnetic field vectors, 
and show each in Figure [7] 


1. We start with that subset of the vectors from Figure[5] which 
satisfy the alignment criterion A 9b < 20°, indicating the location 
of the emissive and absorptive media that we postulate are threaded 
by the same, or nearly the same, magnetic field. This first subset 
is shown in Figure |7Ji, and reveals two features: (i) pixels at the 
extreme NW and SE comers of the CanariCam frame where the 
field seems roughly aligned with the outflow axis from |De Pree et| 
[al~l4 1 994) ; and (ii) mid-frame pixels that seem roughly aligned with 
the orientation of the molecular torus from Howard et al. (1997). 

There is also a narrow trough of low polarisation P close to the 
outflow centre (shown as a long green line in Fig.[7]i), which may 
be due to either: the magnetic field orientation being close to the 
line of sight near this location; an intrinsically weak magnetic field 
here; a field that has structure (whether ordered or chaotic) within 
a resolution element or along the line-of-sight; or a field that is 
still strong, but for some other reason the alignment is perturbed or 
destroyed. The polarisation measurements alone cannot distinguish 
among these possibilities. 


2. We next consider the magnetic field orientation for the 
unaligned pixels (i.e., those pixels failing the alignment criterion 
above; open circles in Fig.|6j. The second subset of vectors, shown 
in Figure [Tjj, reveals the magnetic field orientation due only to the 
absorptive component for these unaligned pixels. This absorptive 
magnetic field is also roughly aligned with the mid-frame field in 
Figure |7|t, and roughly perpendicular to the outflow axis ( |De Preej 
let al.|1994| ). 

Taken together, this consistency in orientation between the 
mid-frame magnetic fields in Figures [7Jt & b strongly suggests a 
generally toroidally-sweeping field across the near inside surface 
of the NIR cavity imaged by Hofmann et al.||2004| , and is consis¬ 
tent with a toroidal geometry of similar orientation to the HCO + 
disk of |Howard et al.|fT997| l. 


3. In contrast, the vectors representing the emissive magnetic 
field (Fig.JTJ:), drawn from the same unaligned pixels as the ab¬ 
sorptive vectors (Fig.JTjt), are very different in orientation from this 
last subset. Figure [7J: shows that the magnetic field in this region 
is again roughly aligned along the outflow axis, perhaps tracing 
the field being carried away by the ionised outflow ( |De Pree et al.| 
11994( t. In fact, the emissive magnetic field is aligned almost N-S 
in the inner regions, before bending to be NW-SE to the north and 
south. This is very reminiscent of the morphology of the ionised 
outflow in Figure 2 ofjD^ Pree et al.]|1994| , which also looks N-S in 
the inner contours and then NW-SE in the outer contours. Also, the 
emissive magnetic field is aligned almost precisely with the outflow 
cavity wall delineated by the X structure (or inverted-V structure) 
in [Hofmann et al. | (|2004| see also Fig.[8] next). 


These detailed features are broadly consistent with the orien¬ 


tation of the larger-scale (5 or 0.2pc) components (Smith et al. 


|2000| ), not only in K3-50, but also in a number of other sources 
according to the analysis of |Wright|p007[ >. For example, the rel¬ 
atively large inclinations between the magnetic fields in the emit¬ 
ting and absorbing regions of large parts of K3-50, and thus prob¬ 
ably between the outflow and torus (see § §4.2. HE .2.2) , possibly 
explains the non-detection of polarisation at 100 pm by Novak et 
al. (1989). In their relatively large beam, and with the polarisation 
effectively being integrated along the line of sight, the two field 
components would cancel each other. This suggests an intimate 
connection between the magnetic fields in K3-50 and the structures 
revealed at other wavelengths. We explore these relationships quan¬ 
titatively next. 
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Figure 5. Overlay of the two components from Fig.[3] again computed for 
each 3x3 binned pixel. This sampling scale (0 / /24x0 / /24) corresponds 
approximately to one resolution element (~0"3 including a modest see¬ 
ing effect). Here the orientation of the emission vectors (blue) is rotated 
by 90° compared to Fig.|3] the absorption vectors (red) are unrotated from 
Fig|3] The underlying cyan contours are the same as the blue 10.3 fim con¬ 
tours in Fig.[3] The rotation of the P em vectors means that both vector sets 
displayed here should represent the projected orientation of the embedding 
magnetic field in each of the absorbing or emitting layers, along the line of 
sight in each (binned) pixel. 


4 DISCUSSION 
4.1 Overall Geometry 

Upon comparing Figure s|T]j7] some general trends are immediately 
obvious. One is that the magnetic field structures are embedded 
in the diffuse emission. The emission peaks do not seem to play 
an important role in distorting the field, although the western I- 
peak (cyan in Fig.|TJ seems to have relatively low polarisation, and 
so may be partially responsible (because of beam dilution) for the 
trough of low polarisation mentioned in j ]3.2| The eastern compact 
/-component (red in Fig.|TJ seems to have a fairly constant angle 
of polarisation in all 3 filters (Fig.[2]>, and is also brightest in I in 
the 12.5/im filter. This suggests that it is a cold object(s) with sig¬ 
nificant silicate absorption. 

The other striking thing is that the vectors in Figure [Tji, where 
the emissive and absorptive components are coaligned, are confined 



0 50 100 150 

Absorptive PA (degrees) 


Figure 6. Pixel-by-pixel plot of position angle differences A 0 = 0 a b s — 0 em 
between the emissive and absorptive components of the silicate feature vs. 
position angles 0 a b s from the absorptive component. These are the same 
data as shown in Fig .[5] except that the emissive component PAs are not 
rotated by 90° in this plot. The error bar shows the typical uncertainties. 
For reference, a horizontal black line is drawn at A0 = 90°, where the 
inferred magnetic field orientations from the two components are aligned 
with each other. One can clearly see that these orientations cluster into two 
distinct groups. Defining A0g = |A0 — 90° |, we draw horizontal dotted 
lines to indicate this separation, and use filled symbols where A 0b < 20°, 
while open symbols have A 0b > 20°. In other words, filled symbols have 
magnetic field orientations in the absorptive and emissive layers aligned 
to within 20° (hereafter the “aligned” points), while open symbols have 
magnetic fields that are more poorly aligned than this (hereafter the “un¬ 
aligned” points). 


to the outer regions of the map, and generally avoid the bright struc¬ 
tures. This presumably means that the dust in these regions is rela¬ 
tively cold and again the polarisation is probably being traced in the 
diffuse material. This suggests that it is necessary to measure po¬ 
larisation in the faint diffuse structures rather than the bright peaks 
to elucidate the dominant field morphology. 

Building on the previous work of De Pree et al.| ( |1994[ >, 
|Floward et al.| \\991) , |Hofmann et al.| ( 2004| ), in Figure [8] we 
schematically show, as a set of overlays on the |Hofmann et al.| 
< j2004j > NIR image, the alignment of our polarisation maps with 
several of the salient features of K3-50. Based on these compar¬ 
isons, we have constructed a toy model for the arrangement of 
features in this source, and especially showing the relationship of 
the derived magnetic field configurations to the observed phenom¬ 
ena; this is shown in Figure [9] We note especially the alignment of 
the two main magnetic structures with other features thought to be 
common in star forming regions (namely, disks and outflows). The 
toroidal magnetic field associated with polarised silicate absorption 
is aligned with a structure that is likely to be dominated by a large 
scale rotation: therefore, it suggests that the magnetic field here is 
also rotation-dominated (by this we mean that the rotation of the 
gas in the torus probably drags the field mostly along with it, and 
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Figure 7. Maps of magnetic field orientations. In all 3 panels, the dotted black contours here are the same as the blue/cyan Si4 polarised intensity contours 
in Figs.|2][3] and [5] (a) Magnetic field orientations from that subset of pixels in Figs .[5] and [6] where the field giving rise to both the polarised absorption 
and emission is aligned, i.e., A9g<20°. The location of a narrow trough of low polarisation P close to the outflow centre is shown by a green line, (b) 
Absorptive-only magnetic field orientations from pixels with unaligned fields, (c) Emissive-only magnetic field orientations from pixels with unaligned 
fields. 


stretches the field lines enough to impart some mean field proper¬ 
ties). Similarly, the axial polarised-emission-related magnetic field 
is associated with a strong outflow of material from the central 
source(s) of this region, suggesting that its geometry too is deter¬ 
mined by the outflow. However, the opposite may also be true, i.e., 
that the magnetic field is actually dominant in both domains, and 
controlling the motions seen by |De Pree et al.| ( [T9 94) and Howard] 
et al. (19971. Clearly, we need a way of discriminating between 
these scenarios: whether the magnetic fields in K3-50 are the driv¬ 
ing agent for the structures and motions we see, or whether the 
fields are merely “passengers” as other processes control events. 


4.2 Physics of Observed Structures 

In order to determine the relative importance of the magnetic field 
in K3-50, we first need to ascertain the role played by other physical 
effects, such as gravity, rotation, thermal pressure, and turbulence. 
|De Pree et al.|(1994^ and |Howard et al.] ( |1997| > provide excellent 
starting points for estimates of cloud mass, gravity, rotation, and 
other bulk parameters. 


4.2.1 The Molecular Torus 

We consider first the molecular torus of |Howard et al.| ( fl997| ). Their 
self-consistent model for the torus’ mass and rotation (which is ap¬ 
parently close to Keplerian) based on their HCO + data gave a best 
fit relation of 

Vrot (77) = 2.0 km s -1 (77/0.10 pc) 05 , (1) 

from which one can infer a gravitational mass contained within ra¬ 
dius R of 

M rot (< R) = V 2 R/G = 93M 0 (77/0.10pc) 2 . (2) 

For example, this corresponds to a mass of ~450Mq within the 
nominal 0.22 pc radius of the inner torus, and a total cloud mass 
~5200Mq within 0.75 pc[] Because their mass estimates match 
reasonably well the mass profile derived from the HCO + integrated 
intensity, jHoward et al.jfl997| concluded that the molecular torus’ 
motions are dominated by rotation. From the degree of mismatch 

1 This mass is --double the value quoted by Howard et al. i 1997) for their 
integration of the same model. The different mass estimates may be recon¬ 
ciled if their assumed HCO - abundance is made half as large. 
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Figure 8. Main features at the centre of K3-50. The underlying 2 fim image is from |Hofmann et al.|j2004} , and shows the possible central ~0.1 pc wide 
cluster of IR sources inside the southern portion of a biconical cavity (suggested by the X-shaped boundary around the IR peak) of size ~0.2 pc. The green 
line indicates the approximate direction of the outflow axis, which we take as N158?5E jDe Pree et al.|1994[|Howard et al.jl997| . The black ellipse shows 
the approximate size (0.44 pc) and orientation (minor axis at —21? 5 from N) of the HCO r toroid, also from Howard et al. (1997). The overlays of coloured 
polarisation vectors are the same as in Fig.[7] At a distance of 8.7 kpc, the scale is 20 pixels = l'/6 = 14,000 AU = 0.067 pc. 


between [Howard et al.| ( 1997 )’s model and pure rotation (mostly in 
the inner torus, at radii 0.1-0.22 pc), we estimate that other effects 
(i.e., turbulence, magnetic fields, thermal pressure) could contribute 
as much as ~20% to the observed HCO + kinematics, if one as¬ 
sumes that the |Howard et al.| i j 1997} model is accurate. We explore 
these contributions next. 

|Howard et al.| ( | 1 99~7[ estimated the gas temperature in the torus 
to be 45 K. We compare the contribution of thermal energy in the 
molecular gas to the above estimates, via the Virial Theorem: 


Af th 


okTR 

mn 2 G 



(3) 


Again, at a radius of 0.22 pc, this gas temperature would provide a 
thermal energy content sufficient to balance the gravitational self¬ 
energy of a ~47 Mq gas cloud. Since this is much less (~25% at 
0.1 pc, ~10% at 0.22pc, ~3% at 0.75pc) than the mass inferred 
from rotational-gravity balance at this radius, it indicates that ther¬ 
mal energy in the gas contributes a relatively small fraction to the 
energy content (or dynamics) of the cloud in the torus, and that this 
fraction drops off to even smaller contributions at larger radii, as¬ 
suming the temperature is constant. If the gas temperature actually 


drops with radius, as seems likely, then the thermal contribution to 
the energy content is even less. 

Similarly, the Virial term giving the mass in the torus that can 
be supported against gravity by magnetic pressure is 

(4) 

At this point, we have no estimate of the magnetic field strength^] 
B in K3-50. Nevertheless, if we assume that the magnetic field 
strength is related to the gas density according to the empirical for¬ 
mula] 

2 We use SI units throughout this paper, including tesla instead of gauss for 
the magnetic field strength B. While 1 nT = 10 /iG simply converts the field 
strength, many astronomical formulae are based on cgs units, which do not 
explicitly account for the vacuum permeability fio . We have made all these 
adjustments in the formulae an d text. 

3 Although eq. (5) from the |Crutcher|{2012} compilation of clouds may not 
necessarily apply to selections of positions within a single cloud, our intent 
here is merely to demonstrate how one scaling argument might work in this 
case. We show in j ]4.2.3| that the estimates made here may actually be quite 
reasonable. 
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Figure 9. Schematic di agram of essen tial features of K3-50, shown as idealised geometries with approximate sizes, shapes, and orientations. In grey we show 
the HCO“*~ toroid from [Howard et af]| 1997f , with the major axis size as indicated and the toroid’s rotation from the HCO+ data also shown by the red and 
blue arrows at the toroid’s ansae. The inner part of the IR cluster of sources from |Hofmann et ah] j2004) is shown in gold, while the fainter IR cavity around 
this appears in orange. Using our CanariCam data, we discern from the polarised silicate absorption a toroidally-configured foreground magnetic field (green 
arrows) that seems to share the toroid’s geometry and perhaps also its rotation. In contrast, the magnetic field (blue arrows) giving rise to the polarised silicate 
emission seems to be aligned with the blueshifted motions of the near side of the bipolar ionised outflow (light blue wedge) of[5F Pree et ah](1994) . The 
redshifted, far side of this outflow is partially obscured by strong IR extinction (light red wedge), presumably due to the material in the front side of the HCO - * - 
toroid. The axis of symmetry of this geometry on the sky is rotated -21? 5 from N, as shown to the right of the diagram. On the far left of the figure, and keyed 
to the diagram, we also show the vectors of the observed magnetic field associated with silicate absorption (i.e., a combination of FigsJTJi and b); on the far 
right we show the magnetic field associated with silicate emission (Fig.JTJ;). 


Btot = 49 nT ( n / io 11 !!!- 3 ) 0 ' 65 ( 5 ) 

from Zeeman observations of similar clouds (see Fig. 6 of|Crutcher| 
1 201 2) , and scale it to the |Floward et al.|fl997) density model at the 
inner radius of the torus 

n = 5.0 x 10 n m“ 3 (fi/O.lOpc)" 1 ' 5 , (6) 


we obtain F?tot = 140nT at the inner edge of the torus (R = 
0.10pc), dropping to 3.8 nT at the outer edge of the cloud (R = 
0.75 pc). Combining eqs. (4-6), 

M mag = 270M 0 (i?./0.10pc) 1 025 . (7) 


Comparing this with eq. (2), we see that magnetic pressure at this 
level could provide ~300%, 133%, or 40% of the support required 
to balance gravity at 0.1 pc, 0.22 pc, or 0.75 pc, respectively. As 
noted above, while estimates based on eq. (7) do not arise from ac¬ 
tual field strength measurements in K3-50, they do make reasonable 
the hypothesis that magnetic fields may contribute significant sup¬ 
port against gravity at the inner edge of the torus. At the same time, 
eq. (7) suggests that magnetic support is less likely to be significant 
at the outer edge of the cloud. 

Another way to write eq. (4) is through the mass-to-flux ratio 
l |Crutcher et al. 12004} , 


(M/<f>) obs 

(M/$)crit 


0.064 


A^H 2 /10 24 m- 2 

BTOf/nT 


( 8 ) 


[Here we use = 2.35 for the mean molecular weight, from a 
slightly more accurate value of 9% for the He fraction by number in 


the gas, based on recombination-line studies of HII regions across 
the Galaxy ( [Shaver et al.|1983] >, instead of Crutcher et al’s fi = 2.8, 
which appears to be a typographical error, for 10% He. This makes 
the numerical coefficient in eq. (8) ~16% smaller than injCrutcherj 
|etaL| l |2004) .] 

In eq. (8), A>1 means that gravity dominates over magnetic 
pressure support, and the cloud is considered “supercritical.” Using 
the above estimates for B, n, and R at the inner and outer edges of 
the molecular torus, we find that Ai nn er~0.68 while A ou ter~140. 
Again, we see that gravity may be slightly “subcritical” (i.e., mag¬ 
netic support dominant) near the inner edge of the torus, but that 
near the outer edge the situation is strongly supercritical (i.e., grav¬ 
ity dominates). 

We caution that the above estimates assume that the B-n scal¬ 
ing relation in eq. (5) applies in the K3-50 torus. We should com¬ 
pare actual magnetic field strength measurements to gauge if these 
estimates are realistic. The only existing Zeeman measurements in 
K3-50 are for OH maser spots located near the NE comer of the 
UCHII region and IR cavity (e.g., [Fish et al.[2005[ . which is also 
near the interior NE ansa of the molecular torus, and give magnetic 
field strengths ranging over 300-700 nT. While indicating strong 
fields, the very restricted distribution of these maser spots lim¬ 
its their relevance to the thermal (whether molecular or ionised) 
gas. If K3-50 becomes a future subject of a Zeeman experiment at 
~arcsecond or better resolution using the dense molecular tracer 
CN, then we would be in a better position to examine the physical 
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roles of all these processes. However, we explore below some direct Thus, the expectation is that the kinetic energy of the ionised 

magnetic field estimates from our polarisation measurements. outflow dominates all other physical effects therein, to be consid¬ 

ered below. 


4.2.2 The Ionised Outflow 


The situation in the ionised outflow is very different, since the gas is 
much hotter and there is an entirely different bulk motion in the gas. 
|De Pree et al.|(T994| presented 14.7 GHz continuum and recombi¬ 
nation line observations of a roughly cylindrically-shaped, bipolar 
ionised outflow extending 15” north and south from an ultracom¬ 
pact HII region, and derived an outflow speed of V = 35 km s -1 . 

Based on their data and model for the HII region, we infer a 
median electron density n e in the outflow area of ~3 x 10 9 m -3 , al¬ 
though this drops quickly as one goes away from the UCHII peak, 
where n e ~ 8xl0 11 m -3 , down the l = 0.63pc-long and 2r = 
0.21 pc-wide outflow channels (measured from the source centre). 
Thus, a rough estimate for the median kinetic energy density in the 
ionised outflow is 

/ \MV 2 \ , a , 

E kin = 2 ( )=P eV = 8 ' 6 x 10 J m ’ (9) 


where we have used the median n e quoted above to calculate the 
equivalent mass density p e = p rriH ti e and taken the mean mass per 
particle p = 1.4 (now in the ionised, not molecular, gas). This esti¬ 
mate for -Ekin scales with n e , so even though the bulk motion near 
the UCHII peak seems less than the average speed qu oted above 
(~10 kms -1 from Figs. 5 & 6 of De Pree et al. 1994 i, -E k i n may 
be several times larger than the value in eq. (9) at the UCHII peak. 
What is certain is that this kinetic energy density is higher closer to 
the outflow centre, and lower at the far ends of the outflow. 

We may similarly infer a thermal energy density in the HII 
region, using 

E th = ^p e o 2 = n e kT e = 3.3 x Id -10 J m -3 , (10) 


where a is the one-dimensional velocity dispersion in the HII 
plasma with a measured electron temperature T e = 7900 K from 
|De Pree et al.HISgg) . Again, we have used our median value for 
n e as above as a rough estimate, but we can see that this is smaller 
than the kinetic energy density in the bulk outflow, by a factor of 
~30. Since the ratio of energy densities EkinAEth = § {V/o) 2 is 
independent of the particle density, it is clear that the kinetic en¬ 
ergy everywhere dominates the thermal energy in the HII region. 
In other words, since the outflow is supersonic (M ~ 3), it should 
not be surprising that the dynamical pressure exceeds the thermal 
pressure. 

These calculations and estimates would ideally be compared 
to measures of the magnetic energy density in the gas (whether 
molecular or ionised), given by B 2 /2po (where po is the perme¬ 
ability of free space), e.g., through Zeeman polarisation measure¬ 
ments of Bio S at cm or mm wavelengths. We describe below some 
magnetic field strength measurements for K3-50, but for now, we 
again estimate what magnetic field strengths would be energetically 
significant in the above settings. From eq. (9), we find that any 
magnetic field in the outflow would need to rise to a level where 
B 2 /2po ~ 10 -8 J m -3 , or B ~ 150nT, in order to be energeti¬ 
cally significant. Such a magnetic field in the thermal ionised out- 
flow/UCHII would be unusual (i.e., discounting the several x 100 
nT fields in the OH maser spots): typical measures in HII regions 
from (e.g.) OH Zeeman absorption observations in foreground cold 
molecular gas are 0.3-3 nT (e.g., 3.8 nT in NGC 2024;'Barnes et ah 
[T989l|Crutcher|2012) . 


4.2.3 Chandrasekhar-Fermi Methods 

Despite the lack of Zeeman magnetic field strength measurements 
for K3-50, we do have a rather unique set of MIR polarisation posi¬ 
tion angles (f? a bs and 9 C m ), which can support some detailed anal¬ 
ysis. 

The standard Chandrasekhar-Fermi analysis l jChandrasekhar| 
|& Fermi 1 1953 [ i connects the dispersions in polarisation angle, 89 = 
s (radian), and velocity, 8V (cm s -1 ), to the plane-of-sky magnetic 
field strength, E pos (G). According to Crutcher et al7| i |20()4} , in 
these units the relation is 

Epos = Q y/Ur p SV/s ; 

making the SI conversions, we have 

Epos = 0.85 pT y/n AT/s , (11) 

where Q = 0.5 is a numerical factor to correct for various smooth¬ 
ing effects (e.g., see |Ostriker et al.|2001}, p = p man as before (n 
in m -3 ), the velocity FWHM AV' = \/81n2 8V (kms -1 ), and s 
is now evaluated in degrees. [Again, we use a smaller value of p 
appropriate for a slightly smaller He fraction.] If we take 5 V ~ 
3 kms -1 as a median value in the HCO + torus from lHoward et all 
< | 1997) > and n from eq. (6), then we will be able to relate s in our 
polarisation maps to the field strength E. We now examine some 
statistics of the dispersion before returning to this standard analy¬ 
sis. 

One approach to evaluating the behaviour of s is that ofjMyersj 
[& Goodman| ( |1991| l. In their language, the goal is to identify a “cor¬ 
relation length” in the implied magnetic field orientation, within 
which the magnetic field directions are correlated and aligned with 
each other, and outside of which they are not. 

Using the formalism of |Myers & Goodman| i j 199~i~j l, we fit the 
distributions of polarisation position angle 9 with a simple gaus- 
sian e b' to obtain a best-fit value for the dispersion s in 
9b (measured in radians). This is a simplified version of Myers 
& Goodman ( [ 199 l| >’s analysis, since they showed that this ap¬ 
proach gives very reliable results even for their comprehensive data 
(i.e., ^hundreds of stellar polarisation measurements) on the Tau¬ 
rus molecular clouds. We have a smaller data set of 9b, so will not 
need the full |Myers & Goodman| l [T99T[ l treatment. 

To evaluate s in our data, we first construct histograms of the 
position angles from various subsets of the data as displayed in Fig¬ 
ures |5jj7] The first is # a b s for all the vectors that can be associated 
with the silicate absorption in the torus as in Figures [TJi and b, and 
is given in Figure [TO] (The analysis for the emissive component in 
Figure[TT]is slightly more complex, so we discuss it after first un¬ 
derstanding the absorptive analysis.) We consider only the polarisa¬ 
tion angle distributions within a sub-window of our data (illustrated 
in Figs.[T()|and[TT|, where we can highlight the differences between 
the absorptive and emissive components. Note that all mean quan¬ 
tities among these pixels that we calculate below are unweighted in 
our analysis, since the main sub-window that we examine here was 
also chosen to cover only high S/N pixels. 

In the absorptive case, the distribution of 0 a b s is well- 
approximated by a gaussian, with a mean position angle (indicat¬ 
ing the large-scale orientation of the magnetic field, in this case 
apparently aligned with the HCO + torus) and dispersion (which is 


© 2015 RAS, MNRAS 000,[Tp?? 








































12 Barnes et al. 
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Emissive B-field 
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P.A. (radian) 

Figure 10. (top) Embedding magnetic field orientation for P a bsi these are 
the same red vectors indicating $b a b s for the absorptive component as 
in Fig.Bl (bottom) Histogram of #B,abs at all pixels within the indicated 
2'/5x2 7 ' box in the top panel. Such histograms can also be constructed for 
any smaller area within the box shown. See text for further discussion. 



P.A. (radian) 

Figure 11. Like Fig. |10| but for the emissive component of the magnetic 
field, (top) Embedding magnetic field orientation for P em ; these are the 
same blue vectors indicating $B,em for the emissive component as in Fig.[5] 
(bottom) Two overlaid histograms of $b em at all pixels within the indi¬ 
cated boxes in the top panel. The black histogram is for all pixels within the 
rectangular black box, while the green-hatched histogram is for all pixels 
within the smaller, irregular green box that includes the higher-polarisation 
pixels on the western side of the black box. Such histograms can also be 
constructed for any smaller areas within the boxes shown. See text for fur¬ 
ther discussion. 
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a combination of the small measurement uncertainties, and the in¬ 
trinsic degree of disorder in the large-scale magnetic field). We first 
focus on the dispersion sg a in such # a b s distributions. For example, 
in the histogram of Figure [lO] the dispersion sg a is approximately 
0.35 radian. 

We can also construct such histograms for subsets of this area, 
and then compute a dispersion in the # a bs distribution for each such 
smaller box. The smaller these boxes are, the more choice we have 
of where to fit them into the large box. Therefore, we can obtain 
a mean ± standard deviation in the dispersion <sg a (A)> for all 
boxes of a given area A. Finally, we plot these results for each box 
size, ranging from a minimal useful size of 3 x 3 pixels (correspond¬ 
ing to 9x9 pixels in unbinned maps), to the full box size shown in 
Figure [lO] obtaining the run of mean ± SD dispersions <sg a > as 
a function of the chosen box size (A 1 ' 2 ). 

This result for the absorptive component is plotted in red in 
Figure [12] Flere we have statistically subtracted, in quadrature, the 
contribution of the 6 measurement uncertainties to the measured 
average dispersions <se a (A)>. For the magnetic field associated 
with absorption, s is quite small at all box sizes, but slowly rises 
with A. We attribute this to the (fairly self-evident) high degree of 
order in the absorptive polarisation vectors in Figure [T]}; in other 
words, this polarisation map shows essentially a single “correlation 
length” across the imaged area. This makes sense in our interpreta¬ 
tion of this magnetic field being associated with the torus, since the 
torus is ~3 x wider than our map. 

We can repeat this exercise for the 9 em distribution, except that 
we must first caution that the original CF method §Chandrasekhar| 
|& Fermi 1 1 953| i was not developed for magnetised gases in super¬ 
sonic motion. Nevertheless, numerical simulations of supersonic 
(M = 5-9) MF1D turbulence i jOstriker et al.|200f| show that it may 
still be reasonable to employ the CF approach to estimating mag¬ 
netic field strengths in the outflow, where M ~ 3 (see §4. 2. 2} . 

In this case, the 9 em position angles are clearly not distributed 
around a single mean value. To see this, note in Figure[TT]the dis¬ 
tinct change in <9 em > between the northern and southern portions 
of the green box. Note also that 0 Gm for pixels outside the green 
box, but inside the black box, are oriented more randomly than the 
pixels inside the green box. Because of this behaviour, we perform 
the above analysis for both samples of pixels (i.e., those from both 
the green and black boxes). 

The black histogram in Figure [TT] shows a much more widely 
distributed # em due to the more variable orientations evident in the 
emissive polarisation map. It turns out this is true at all box sizes 
within the black box, so that the plot (shown in blue) of <sg e > 
vs A 1 ' 2 in Figure 112| shows larger values than <se a > at all box 
sizes. The trend with A is also up, but seems to reach a plateau of 
Smax ~ 0.55 at a box size ~0.025 pc = 5,000 AU ~ 0'.'6. 

Now, when we consider the green-box pixels only, we obtain 
an intermediate result. In Figure [T2] for smaller box sizes within 
the green area (plotted in green), we see relatively small <sg e >, 
around the same minimum value of 0.1 radian. In other words, 
on small scales the magnetic field associated with the emission is 
about as well-ordered as the overall magnetic field associated with 
the absorption in the torus. Then as the box size increases, the trend 
of green points rises rapidly to the same plateau in s as the blue 
points. Again, this is clearly due to the fact that at the larger box 
sizes in either the green or black areas of Figure [IT] we are includ¬ 
ing pixels from multiple “correlation lengths” of the magnetic field, 
so the distributions of 6 represented by s must get wider. 

As explained by |Myers & Goodman] <| 1991), where the CF 



Area 1 ' 2 (pc) 


Figure 12. Dispersion of polarisation position angles 6 as a function of 
box size within the areas shown in Figs ) 10| and |11| See text for further 
discussion. 


method applies, these dispersions are related to the ratio of the dis¬ 
ordered vs. ordered magnetic field strengths via 


N 1 / 2 B ± ’ 


( 12 ) 


where N is the number of magnetic field correlation lengths in the 
line of sight, B±_ is the strength of the ordered component of the 
magnetic field, projected into the plane of the sky (i.e., the ordered 
component of f? poa ), and cjb is the dispersion in the strength of 
the random component of B pos . For now, we estimate N from the 
behaviour of s, and so constrain somewhat the ratio <jb!B±. 

From our CanariCam data, we see that typical values for s lie 
in the range ~0.1-0.6, where the lower value is clearly associated 
with single correlation lengths. In those areas of our K3-50 data, 
we find that <jb/B± ~ 0.1-0.3, i.e., the strength of the disordered 
magnetic field is only 10-30% of the strength of the ordered mag¬ 
netic field. If our interpretations are correct so far, this conclusion 
would therefore seem to hold for the whole of the molecular torus, 
and for small areas of size ~ 0.01 pc = 2,000 AU in the ionised out¬ 
flow. This further suggests that, if the magnetic field in these areas 
is dynamically important, it would be the ordered component that 
is most likely to be relevant, rather than the disordered field. But 
this may still be significant: we saw in §4.2.1| (e.g., by comparing 
eqs. 2 and 7) that the relatively ordered magnetic field configura¬ 
tion in the toms implied by our data may be dynamically important 
in the innermost parts (<0.22 pc) of this structure. Actual CF field 
strength estimates in the molecular torus are discussed next. 

On the other hand, the importance of magnetic fields in the 
ionised outflow appears less favourable. Over larger portions of the 
outflow, we see evidence for at least 3 correlation lengths in the 
imaged area, so <jb/B± ~ 0.6\/3 ~ 1. Over smaller portions of the 
outflow, ctbIB± ~ 0.1-0.3. Given that any magnetic field would 
have to be unusually strong in the ionised gas to strongly influence 
the dynamics ( §4.2.2} , it suggests that neither a disordered nor an 
ordered field is likely to achieve this in the outflow, whether on 
smaller or larger scales. But as we shall see, the CF analysis may 
give a surprising result here. 

Returning now to the standard CF analysis described at the 
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beginning of this section, we can evaluate eq. (11) further with these 
values for s. Combining with eq. (6) we obtain for the molecular 
torus 

Bpoa = (710 - 140) nT (.R/O.lOpc) -0 ' 75 (13) 

for areas with s = 6°-30°, respectively. What this means is that, 
towards the inner part (R = 0.1 pc) of the torus (the northern part 
of the box in Fig.|10|>, s ~ 0.3 rad (17°), implying that the plane- 
of-sky magnetic field strength B poa ~ 250 nT, which actually ex¬ 
ceeds the critical value needed to resist gravity ( §4.2.1} by a fac¬ 
tor of about 1.8 (A ~ 0.56). Towards the outer part of the torus 
(R = 0.75 pc), the density drops off, whence B pos ~ 55 nT. This 
is large, but because of the greater mass contained within the larger 
radius, A ~ 15, which is still strongly supercritical (meaning grav¬ 
ity dominates). These measures for A are very similar to the ones 
we estimated in §4.2.1| and suggest that the general picture given 
there is qualitatively correct. Therefore, if the CF analysis presented 
here is valid, the magnetic field at the inner edge of K3-50’s molec¬ 
ular torus is about twice the strength suggested by the B-n relation 
eq. (5), placing it somewhat above the mean relation in Crutcher’s 
(2012) Figure 6. 

For the ionised outflow, we can perform a similar calcula¬ 
tion, but in this case the deduced magnetic field strength will be 
much smaller than in the molecular torus. This is mainly because 
in eq. (11) the particle density is generally much lower in the out¬ 
flow than the torus, and while the velocity dispersion in the flow 
is larger, so is the polarisation angle dispersion. Inserting the num¬ 
bers from the discussion in §4.2.2|into eq. (11), we obtain B pos ~ 
47 nT, which is surprisingly large given the discussion in §4.2.2| 
but less than the value (150nT) required for magnetic dominance. 
Flowever, given the inherent uncertainties involved in the use of CF 
methods, it is possible that in some parts of the ionised outflow 
(e.g., the green box in Fig. m the dispersion s may be low enough 
to indicate possible magnetic dominance in such areas. 

In summary, although there are substantial uncertainties in the 
application of CF methods to estimating magnetic field strengths 
(as embodied by the wide range of Q values having been used in 
various studies), the weight of evidence suggests a consistent pic¬ 
ture of K3-50: the magnetic field is probably strong enough to sup¬ 
port the cloud against gravity in the inner parts of the molecular 
torus (R < 0.2 pc), but is not strong enough to resist gravity out¬ 
side this area. While the magnetic field in the outflow also seems 
somewhat too weak to dominate the overall energetics, it may nev¬ 
ertheless be important in some locations. 


4.3 Comparison with Models 

Theoretical simulations of star formation have been a vigorous area 
of activity for many years. Recent examples show that, in order 
to understand the magnetic field configurations inferred from po¬ 
larisation measurements of massive protostellar objects, it may be 
necessary to model fully three-dimensional geometries that include 
misaligned rotation and magnetic axes viewed from arbitrary per¬ 
spectives (e.g., |Machida et al.|2006[|kataoka et al.|2012||Shinnaga| 
|et al,|201~2| >. Our results for K3-50 can be compared with such mod¬ 
els on the physical scales represented in our maps, namely ~ 1,000- 
20,000 AU. 

That the magnetic field associated with the absorptive com¬ 
ponent of the 10/-im silicate feature seems to be well-aligned with 
the toroidal geometry of '[Howard et al.[fl 997^ 's molecular disk, and 
that the magnetic field strength in the inner 0.2 pc of the torus may 
be dynamically significant, suggests that the magnetic field may 


have played an important role in the evolution of this star-forming 
cloud, or may even have dominated the formation of these struc¬ 
tures. 

On the other hand, [Machida et al.| l j2006[ l showed that one 
could reproduce such features, especially a toroidally-aligned mag¬ 
netic field on scales of ~2000-20,000 AU, only if the cloud is 
rotation-dominated and the magnetic field relatively weak (their 
“C” models). For the strong-field case, their “D” and “F” models 
show that toroidal alignment only occurs on much smaller scales, 
<300AU. In support of this weak-field scenario, our estimates in 
§4.2.1 show that, outside the inner edge of the torus modelled 
by |Howard et al.| §1997} (i.e., on the larger scales that may be 
more closely related to the “initial conditions” for this star-forming 
cloud), the magnetic field is indeed unlikely to be dominant. 

The review of [Wright| ( |2007| > showed that, on the larger 
(pc) scale, magnetic fields are preferentially aligned toroidally in 
(pseudo-)disks, but poloidally along outflow axes, and we see here 
that K3-50 provides a classic example of such alignments. It is 
tempting to go further and attribute these alignments in K3-50 
on the larger scale (>lpc) to the general interstellar magnetic 
field setting the orientation of gas collapse along field lines to 
form the molecular torus; also on the smaller scale (<1 mpc), to 
the magnetically-driven outflow mechanism of |Banerjee & Pudritz| 
( |2007]> and similar models of that ilk. The justification is that, al¬ 
though our CanariCam K3-50 data do not directly probe either of 
these scales, physical laws must still be “obeyed” in connecting 
these scales through the 0.01-0.1 pc range revealed here. 

However, there are some difficulties in making this connec¬ 
tion. For example, optical polarisation data of stars in the direction 
of K3-50 over the relevant distance range §Heiles||2000} suggest 
that the larger-scale magnetic field in the vicinity of K3-50 is not 
so conveniently oriented. Also, as described above, the formation of 
the dense clumps responsible for cluster formation may be magnet¬ 
ically supercritical, at least on the larger scales. Likewise, many of 
the protostellar outflow models tend to be for the lower-mass/single 
star case, and also are typically computed over 10 2-2 yr timescales, 
rather than for a massive cluster over the 10 4 ~ 5 yr timescale of the 
outflow or disk kinematics (§ §4.2. BE .2.2} . Finally, in the ionised 
outflow, it seems somewhat surprising (given past Zeeman mea¬ 
surements in HII regions) that the magnetic field strength may be 
strong enough in some locations to influence the highly supersonic 
motions in this volume ( §4.2.2} . On the other hand, it is certainly 
clear that the magnetic field here is not “primordial” to the cloud’s 
formation, and must be dynamically relevant to the outflow itself. 
This suggests that other magnetic field measurements in HII re¬ 
gions must be attributed with care to the appropriate dynamical 
setting. 

Thus, challenges in understanding these data and interpreta¬ 
tions remain. As we have seen here, new instruments like Canari¬ 
Cam will certainly be pivotal in improving this understanding. 


5 CONCLUSIONS 

Decomposition of magnetic field components using mid-IR imag¬ 
ing polarimetry is a very clever “old” technique, made new and 
more powerful with CanariCam. By carefully considering the po¬ 
larisation information from this new instrument, we can conduct 
pioneering science such as “polarisation tomography” of magnetic 
fields to decipher their three-dimensional geometry in sources, as 
well as constrain magnetic field strengths in these same areas. 

We have presented new imaging polarimetry results using the 
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mid-infrared silicate feature for the luminous UCHII region K3- 
50, which give new constraints on the magnetic field structure, 
strength, and role in this massive star-cluster forming cloud. Our 
main results are: 

1. We find that the respective magnetic field components are 
substantially aligned with both the ~0.4-pc scale molecular torus 
that surrounds the young star cluster, and with the ionised outflow 
that is being driven from the same proto-cluster. 

2. Scaling arguments and comparison with models suggest that 
the magnetic field may be energetically important in the inner parts 
of the molecular torus, but unimportant outside a radius ~0.25 pc, 
where the torus is instead probably rotationally-supported against 
gravity. 

3. In the outflow, the magnetic field is probably not globally 
dominant, but it may be important in some locations, such as in 
confining the outflow along the cavity walls. 

Future field-strength measurements, such as from Zeeman experi¬ 
ments with OH or CN lines, would be very helpful to confirm these 
results. 

Similar observations of other sources with CanariCam will al¬ 
low us to examine what role magnetic fields play in star formation, 
stellar and protoplanetary disks, AGN, and other astrophysical sit¬ 
uations, providing new constraints on theoretical models for such 
objects. 
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